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ABSTRACT 

The epoch of reionization is a milestone of cosmological structure formation, marking the birth of 
the first objects massive enough to yield large numbers of ionizing photons. Ho"wcver, the mechanism 
and timescale of reionization remain largely unkno'wn. Measurements of the CMB Doppler effect from 
ionizing bubbles embedded in large-scale velocity streams - known as the patchy kinetic Sunyaev- 
Zel'dovich (SZ) effect - can be used to constrain the duration of reionization. When combined -with 
large-scale CMB polarization measurements, the evolution of the ionized fraction, Xg, can be inferred. 
Using ne'w multi-frequency data from the South Pole Telescope (SPT), -we sho'w that the ionized 
fraction evolved relatively rapidly. For our basic foreground model, "we find the kinetic SZ pcwer 
sourced by reionization at i — 3000 to be -Dfooo'^'*' — 2.1 /iii'^ at 95% confidence. Using reionization 
simulations, we translate this to a limit on the duration of reionization of Az = Zx^=o.20~Zx^=o.99 4.4 
(95% confidence). We find that this constraint depends on assumptions about the angular correlation 
bet'ween the thermal SZ po'wer and the cosmic infrared background (CIB). Introducing the degree of 

correlation as a free parameter, we find that the limit on kSZ po"wer 'weakens to ^3000'^'*' ^ 4:.9fj,K^, 
implying Az < 7.9 (95% confidence). We combine the SPT constraint on the duration of reionization 
'with the WMAP measurement of the integrated optical depth to probe the cosmic ionization history. 
We find that reionization ended 'with 95% confidence aX z > 7.2 under the assumption of no tSZ- 
CIB correlation, and z > 5.8 "when correlations are allo'wed. Improved constraints from the full SPT 
data set in conjunction 'with upcoming Herschel and Planck data should detect extended reionization 
at > 95% confidence provided Az > 2. These CMB observations complement other observational 
probes of the epoch of reionization such as the redshifted 21 cm line and and narro"w-band surveys for 
Lyman-a emitting galaxies. 

Subject headings: cosmology: theory - intergalactic medium - large scale structure of universe 
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1. INTRODUCTION 

Galaxy-sized dark matter halos first collapse at z > 25. 
The stars and black holes that form in these halos ulti- 
mately ionize and heat the intergalactic medium (IGM). 
However, exactly when and how this process occurred is 
unknown. 

To date, there have been two primary observational 
constraints on the reionization era. First, Lyman-a 
(Lya) forest absorption spectra towards high redshift 
quasars show that the opacity of the IGM to Lya pho - 
tons is rapidly increasing at z > 6 (e.g., Fan et al. 2006b). 
This increase in opacity has been interpreted as evidence 
for an increasing neutral hydrogen fraction. However, the 
interpretation of quasar absorption spectra is hampered 
by the large cross section for Lya absorption, which can 
lead to complete absorption even if the hydrogen ioniza- 
tion fraction is high. At the very least, the amount of 
transmission in the Lya forest at z < 6 indicates that 



the bu lk of reionization occurred at higher redshifts ( Mc- 
Greer et al.|l2011i) 

Second, WMAP measurements of the optical depth 
through large-angle cosmic microwave background 
(CMB) polarization anisotropy suggest that the red- 
shift of reionizatio n assuming an instant aneous process 
is z ~ 10.6 ± 1.4 dKomatsu et al.||201lD . CMB polar- 
ization measuremen ts otter an inte gral c onstraint on th e 
reionization history ( Kogut et al. 2003 ; P age et al.|2007 ) . 
Therefore, the polarization data are fully consistent with 
either instantaneous or extended reionization scenarios. 

Together these observations have been interpreted by 
some as favoring a prolonged reionization epoch ending 
at z ~ 6 (Miralda-Escude 2003; Gen 2003; Fan et al. 



2006b| [Bolton fc Haehnelt,|2007b, ,Wyitlic fc Gen, 2007^ . 
However, the claim that z ~ 6 Ly-a forest spectra probe 
the tail end of reionization is somewhat controversial. 
These data may be consistent with reion ization complet 
ing at a redshift higher than z = 6 (e.g.,|Oh & Furlanetto 



. .10 

2005 , Lidz et al. 2006, , Becker et al.|2007||McGreer ei al 



2011 



et al.i 



Additional model-dependent constraints on reioniza- 
tion, most of them upper limits on the neutral fraction, 
have been derived by other me ans: (1) the size of the 
proximity zone around quasars dWyithe et al. 20051 |Fan 
2006b, bu t see [M csingcr et al. 2004; BoltonT^ 



Haehnelt 2007a| [Lidz et al. 2007, ,MaseUi et al.|,2007] ,, 
(2) the claimed detections or danip ing wing absorption 



from neutra l IGM in quasar spectra (Mesinger & Haiman 
2007) which were used to place lower limits of 



2004 



2070 and 3% on the neu tral fraction at z=6.3 (but see 



Mesinger & Furlanetto 20 08aJ , and more recentl y a lower 
limit of 10% at z = 71 ( [Mortlock et al.|[2011| ; (3) the 
non-detection of intergalactic damping wing abs orption 



in the spect r um of a gamma ray bu rst at z=6.3 (Totani 



et al ||2006[ [McQuinn et al | | 2008| ; and ( 4) the num 
ber density and clustering ot Lya emitters dMalhotra fc| 
Rhoads 2004; Haiman & Gen 2005; Furlanetto et al. 20(W 
Kashikawa et al. .2006,, McQuinn et aL2007a^ ,Mesinger 
Furlanetto| i2008b[ [Uuchi et al.||2010| ^ These constraints 
will be discussed in more detail in ^ and Figure 13 



In this paper, we produce the first constraint on the 
evolution of the ionized fraction using GMB data. We 
use new data from the South Pole Telescope (SPT) to 
constrain the amplitude of the "patchy" kinetic Sunyaev- 



Zel'dovich ( kSZ) signal resultin g from inhom ogeneous 
reionization ([Knox et al. il998; G ruzinov fc Hu|1 998; Sa^ 
tos et al.| 2003VZahn et al. 2005a| IMcQuinn et al.|2005ar 
Uiev et al.| 2007; Fan et al. 2006a|' We interpret this sig- 
nal in the context of reionization models to place the first 
constraints on the duration of the epoc h of reionization 



using GMB data. A companion paper ( Reichardt et al. 
2012 hereafter Rll) describes the SPT data underlying 
this ineasurement in detail. We combine the SPT data 
with WMAP7 measurements of the integrated opacity 
due to reionization. The combination makes it possible 
to constrain the evolution of the ionized fraction, signif- 
icantly limiting the allowed range of reionization scenar- 
ios. 

The structure of the paper is as follows. In Sj2]we out- 
line the physics of the kSZ effect from inhomogeneous 
reionization ("patchy kSZ") and introduce our simula- 
tion scheme. We describe the data in f|3] In i4l we 
set up the model for the GMB and astronomicalfore- 
grounds which we fit to the data. In ^ we review the 
constraint on the integrated optical depth by WMAP7 
and present the SPT constraints on the amplitude of the 
patchy kSZ component. In we present our constraints 
on the evolution of the ionized fraction, and forecast fu- 
ture constraints with the SPT, Planck, and Herschel ex- 
periments. We summarize our results and place them in 
the context of other observations in S|7j 

2. EPOCH OF REIONIZATION AND KINETIC 
SUNYAEV-ZEL'DOVICH EFFECT 

In this section, we review the theory of the kSZ sig- 
nal, and then introduce the reionization simulations and 
discuss the kinetic SZ power used to fit to the data. 

2.1. CMB reionization observables 

Scattering of GMB photons into the line-of-sight in- 
troduces temperature anisotropy through the Doppler 
effect if there are perturbations in the baryon density, 
Ph, or ionization fraction, x^. These scatterings slightly 
change the temperature of the GMB blackbody. The 
GMB Doppler effect is often called the kinetic Sunyaev- 
ZePdovich (kSZ) effect (fSunyaev fc Zerdovich|1970l ' jSmr 
yaev & Zel'dovich 1980) when referring to reionization 
or the nonlinear regime of struct ure form ation, an d the 
Ostriker-Vishniac effect (Ostriker & Vishniac 1986) when 
referring to the linear regime, but we will call it kSZ effect 
throughout. The total contribution to the temperature 
anisotropy from a redshift interval [zi,Z2] is given by: 



kSZ 



COMB 



(n) 



(Tt 

c 



^'^ dx dz 
^ dz{l + z) 



ne(z)e-"(^)n-q, (1) 



where aT is the Thomson scattering cross-section, ^ is 
the comoving line-element and fi is the line of sight unit 
vector. ne(z) is the mean free electron density. 



ne(z) 



Xe{z)pb{z) 

IJ-e-mp 



(2) 



where Xe{z) and Pb(z) are the mean free electron fraction 
and mean baryon density of the universe as functions of 
redshift and /ieWp is the mean mass per electron. We set 
/ie = 1.22, appropriate for singly ionized helium. Due to 
their similar ionization potential, we assume that helium 
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is singly ionized alongside hydrogen, that is, when all 
hydrogen is ionized and all helium is singly ionized = 
1. Note that Xe = 1-07 when helium is fully ionized; our 
model for the post-reionization kinetic SZ signal assumes 
that this occurs at z — i. 

The optical depth, t[z), from the observer to redshift 
z is given by. 



ne{z') dx 



1 



dz'. 



(3) 



The WMAP large-scale polarization signal constrains the 
optical depth to the end of recombination, t(z ~ 1000) — 
0.088 ±0.014. 

Finally, inhomogeneities in the ionization fraction and 
baryon density enter through 



q= (l + 4)(l + 5b)v. 



(4) 



where 5x = Xe/xe — \, 5b = Phi Pb — 1, and v is the bulk 
motion of free electrons with respect to the CMB. 



We discuss the post-reionization kSZ signal in { 4.1.1 
Reionization produces kSZ power primarily through loca^ 
changes in the ionization fraction. To predict the ampli- 
tude of this signal we require a prescription for calcu- 
lating the evolution of the ionization morphology during 
the epoch of reionization. 

2.2. Efficient Monte Carlo reionization simulations 

Our ability to infer the properties of the first ionizing 
sources from observations hinges on the accuracy with 
which the ionization morphology during reionization can 



diative transfer codes (e.g., fGncdin 2000; Sokasian et al. 


2001; Razoumov et al.||2002; 
et al.,|2006 |McQuinn et al., 


Ciardi ct al. 2003; Mcllema 
2007b, .Semehn et al..,2007r 


Trac &: Cen||2007al |Altay et al.||2008l |Aubert & Teyssier 


2008 Finlator et al. 2009 


Petkova & Springel 2009). 



for simulations of reionization are daunting. Simula- 
tions must resolve the low-mass galaxies (down to the 
atomic cooling threshold corresponding to a host-halo 
mass ~ 10* Mq at z ~7-10) that are expected to domi- 
nate the ionizing photon budget. They must also be large 
enough to statistically sample the distribution of ionized 
regions. These regions can span tens of comoving Mpc 



in size towards the end of reionization (F urlanetto fc Oh 
20051 IZahn et all|2005b[ [20071 fMesinger fe Jj'urlanetto 



2007 



Shin et al.|2008p . Recently, some groups have come 



close to achieving this dynamic range in a single simu- 
lation dlliev et al.|l2006: McQuinn ct al. 2007b; Trac & 



Cen"2007a';'yhin et al. 2008; Trac et al.|2008(|Zahn et al. 



2011,; sec Trac & Gncdin 2011 for a recent review) 



Given the large uncertainties in the production and 
escape rates of ionizing photons in high-redshift galax- 
ies (and also in how the photons are absorbed by dense 
systems), a large parameter space must be explored to 
interpret observations. These concerns have prompted 

^ Note however that Lyman hmit systems (LLSs), which can 
dominate th e absorptio n of ionizing photons (see for example the 
appendix of [Furlanetto fc Oh 2005} , are still too small to be re- 



cally included via analytic prescriptions (e.g., 


Furlanetto & Oh 


2005 


IMcQuinn et al.|2007b||Choudhury et al.|2a 


U9||lJrociam et al. 


2011 





several groups to deve lop much more CPU efficient, ap- 
proximate algorithms (|Furlanetto et al. 120041 iZahn et al 



|2005b','200 7ljMesinger fc Pur ianetto'i^OO?^ 'Geil k Wyithe 
,2008 Alv arez et al.|20091|Gho udhury et al. 2009; Thomas 
|et al.|2009 Zahn et al.|2011 ) . The basic idea is that semi- 
anaiytic rnodels can be apiDlied in a Monte Carlo fash- 
ion to large-scale realizations of the density and veloc- 
ity fields. The resulting complex reionization morphol- 
ogy can be compared side-by-side with simulations based 
on the same initial conditions combined with the same 
source and sink prescr iptions. We use a modification of 
the I Zahn et al.] ( 201 1[ ) model to set up our reionization 
parameter space j 

The model associates ionized regions with the ionizing 
sources they contain. The size distribution of ionized re- 
gions is related to the halo mass function through the 
ansatz, Mjon = C-^gai- Mg^i is the mass in collapsed ob- 
jects. C is the efficiency factor for ionization, and can 
be decomposed as C = fcscf*N^/bn~^^. Here, /esc is the 
escape fraction of ionizing photons from the object, 
is the star formation efficiency, N^/b th^ number of ion- 
izing photons produced per baryon converted into stars, 
and riioc is the typical number of times a hydrogen atom 
has recombined. Although the efficiency factor is a rough 
combination of uncertain source properties, it can encap- 
sulate a wide variety of reionization scenarios. 
We calculate M ^^,] according to the extended Press 



Schechter model (Bond et al. 1991). In this model, the 



collapsed fraction (or the traction of baryons that lie in 
galaxies) in a region of size r depends on the mean over- 
density of that region, S{r), as: 



/cou(f, Mmin) = erfc 



6c{z)-Sir) 
v/2[a2(r„,i,)-a2(r)] 



(5) 



Here, Sc{z) is a numerical factor from linear theory equal 
to 1.686 today, and ^^(r) is the linear theory rms fluctu- 
ation smoothed on scale r. rmin is defined as the radius 
that encloses the mass Mmin (at average density p) cor- 
responding to a virial temperature of 10* K, above which 
atomic hydrogen line cooling becomes efficie nt. The red 



shift depen dence of the minimum mass (see Barkana & 
Loeb|[200T ) is given in terms of the virial temperature: 



M, 



10« 



10 



2 X lO^K 1 + z 



3/2 



no A,(z) 



-1/2 



(z) 187r2 

(6) 

with the fitting function Ac(z) = 18tt^ + 82d - 39d^ 
(Bryan & Norman 1998), where d = flm{z) — 1 is eval- 
uated at the collapse redshift. The luminous mass in 
galaxies required to fully ionize all hydrogen atoms is 
inversely proportional to the ionizing efficiency, so ion- 
ization requires 



/coll > C 



(7) 



Hence, we can define a barrier Sx{r, z) which fluctuations 
have to cross for their baryonic content to become ion- 

^ The reason we do not use purely analytic models is that they 
have difficulties d escribing the partial ove rlap stage of merging ion- 
ized regions (e.g., [McQuinn et al.|2005b| . 
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ized: 

S{r) > 5,{r,z) = 5,{z)-V2cric-\C^)[a\n,,,,,)-a\r)Yl 

(8) 

The algorithm then proceeds as follows. For every po- 
sition, the linear matter over-density is calculated within 
a range given by the smoothing kernel (initially set to 
a larg e value comparable to the simulation size) . iZahn 



et al.| (20111 show that good agreement with radiative 



transfer simulations is achieved when the smoothing ker- 
nel is a top hat in harmonic space. The halo collapse frac- 
tion is estimated from this over-density, and translated 
with the efficiency factor Q into an expected number of 
ionizing photons based on the number of collapsed halos 
above the atomic cooling threshold mass. The number 
of ionizing photons is then compared to the number of 
hydrogen atoms at each point. If there are sufficient ion- 
izing photons, a point is labeled ionized. If not, a smaller 
radius is set and the algorithm repeated until the reso- 
lution of the simulation box is reached. 

We optimize the algorithm for mock patchy kSZ spec- 
tra at the SPT angular resolution. First, we simulate 
the density and velocity fields in a rectangular box with 
a high dynamic range, 1^ x 3 Gpc/h with 512^ x 1536 vol- 
ume elements. The simulation box subtends 8.5 degrees 
on the sky at the central z=9 plane, yielding arcminute 
resolution and low sample variance on the scales of in- 
terest. Along the z-axis, the comoving distance is trans- 
lated to redshift. The volume thus extends from 2; ~ 4 
to z ~ 27. The large redshift extent allows us to cap- 
ture most reasonable reionization scenarios. The densi- 
ties, velocities, expansion rate, minimum virial mass for 
star formation, and excursion set barrier (Eq. [8| are ad- 
justed with redshift. The large angular and line-of-sight 
coverage are important because the ionized regions and 
large-scale velocity streams that source the kSZ signal 
appear on scales of tens to hundreds of megaparsecs. 

In addition to the global ionization efficiency param- 
eter, which sets the timing of reionization, it is im- 
portant to allow for an additional degree of freedom to 
capture the influence of physical processes that may im- 
pact the duration^ such as feedback effects and recombi- 
nations. To do this we add a "feedback" parameter, a, 
that makes the effective ionization efficiency an inverse 
power law of the expected ionization fraction in absence 
of feedback, x* , which is a measure of the ionizing pho- 
tons available at a given time|^ 



C = Co 



X* = Co /coll (f = 00, Af„ii„) . 



(9) 



A database of simulated models is set up for C, G 
[10, 1000] and a € [-50, 1]. This leads to - 10000 simu- 
lations that have reionization durations of Az € [0.1, 12], 
and ending redshifts of z^nd G [4, 16]^ Here we have 



defined Az 



eg 



^cnd 1 ^cnd 



Za- 



= 0.99: 



and Zbo 



Zx^=o.20- The choice of a relatively large ionization frac- 
tion of Xe = 0.20 for Zbcg is due to the fact that lower 

^ For a constant ionization efficiency C, (« = 0), the simulation 
predicts a reionization duration of Az ~ 3 — 4, given the definition 
of Az below. 

Note that the Ly-Q forest d ata effectively rule out models that 
extend past Zenj < 5 (see e.g., [Fan et al.j|2006b| |McGr eer et al.| 
|2011| , and we will impose this as a hard prior wfien^deriving reion- 
ization constraints in il6l 



ionized fractions are largely made up of ionized regions 
too small to be probed by the current SPT data, hence 
we do not attempt to constrain the reionization state 
preceding this stage. An example of a 20 Mpc/h deep 
projection of an extended reionization simulation with 

— 12, a — 0.8, and Az ~ 9 is shown in Figure [T] 

For each reionization simulation, we interpolate rays 
from constant comoving positions to constant angular 
positions, produce a map of the kinetic SZ with a size of 
~ 7x 7 degrees, and measure its angular power spectrum. 
Throughout this paper, power spectra are given in terms 
of Di = Ci£{i + l)/27r IfiK'^). For each simulation, we 
also calculate the total integrated optical depth assuming 
helium is singly ionized alongside hydrogen and doubly 
ionized at z=3 (the latter assumption has a negligible 
influence on the results). 

The spatial resolution and box size used lead to con- 
vergence of the patchy kSZ amplitude on the relevant 
angular scale (see the next section). However, the ion- 
ized fraction is underestimated during the early stages of 
reionization, when ionized bubbles smaller than the spa- 
tial resolution of the simulation contribute significantly. 
To correctly predict the duration and integrated opti- 
cal depth we therefore scale the ionized fraction at each 
fraction to the analytical value C/cou('' = oo,M„iin), the 
mean ionization degree of the universe at that redshift. 

In these simulations, we assume a ACDM cosmology 
with nah"^ = 0.111, rJb/i^ = 0.0222, r^A = 0. 736, ns = 
0.96, r 



0.085, and as = 0.8, consistent with Komatsu] 



et al. (2011). We di scuss the cosmological scaling of the 
kSZ predictions in §4.1. 1| and the effec ts of uncertainty 
in the cosmological parameters in { 5.2 



2.3. The kinetic SZ signal from reionization 

The amplitude and shape of the patchy kSZ power 
spectrum depend on the redshift evolution of the ion- 
ized fraction and the reionization morphology. Here we 
discuss the dependence of the kSZ power spectrum on the 
physics of reionization and motivate the template shape 
used to fit to the SPT data in order to constrain the 
ionization history. 

The reionization model outlined in |2.2| has been shown 
to agree well with full radiative transfer simulations on 
scales up to 100 Mpc/h (see e.g., |Zahn et al.|2007[[20TT ). 
As a generic consequence of the increase m source col- 
lapse fraction with time, the size distribution of ionized 
regions and volume fllling factor evolve more slowly at 
the beginning of reionization than at the end. Also, the 
roles of feedback effects and recombinations depend on 
the degree of ionization. If we instead used a model with 
a parametrized yet physically unmotivated analytic form 
of Xe , the relation between the redshift corresponding to 
various ionization fractions Zj^ , the kSZ power, and in- 
tegrated optical depth r would be different, potentially 
biasing constraints on Xe- 

Using radiative transfer simulations, it has been found 
that the size distribution of ionized regions at flxed Xg 
is relatively robust to redshift tra nslations zi — > Z2 (e.g., 
Zahn et al.]|2007[ [McQuinn et al. 2007b ). In addition, 
the angular size corresponding to a given co- moving scale 
varies little across the redshifts of interest. The change in 
the patchy kSZ power spectrum due to a translation from 



7.5 to Zinid — 11.5, where Zmid 



=0.5 



for sim- 



ilar Az is shown in the black solid and dot-dashed lines 
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.pdf 



z=4 z=6 z=10 z=15 z=27 




1.5 Gyr 1.0 Gyr 0.5 Gyr 0.25 Gyr 0.15 Gyr 



Fig. 1. — A 20 Mpc/h deep cut through one of our simulation volumes with dimensions 1x3 Gpc/h, corresponding to an angular extent 
of 8.5 X 8.5 degrees at the central z=9 plane, and a redshift coverage from z ~ 27 to 2: ~ 4. Black means fully neutral and white means 
fully ionized. Note that this is not a result from fitting to data but a random sample from our simulation database. The ionized fraction 
in this particular simulation evolves from zero at z ~ 20 to fully ionized at z=5.6. 



in Figure [2] The shapes of both power spectra are very 
similar, with the main difference being in their normal- 
ization. The amplitude is increased in the Zj^id — 11.5 
model due to the higher mean density of the universe 
at higher redshift, which is counteracted somewhat by 
linear growth of density and velocity fields toward later 
times. 

When moving to very high redshifts, ever rarer sources 
have to yield enough ionizing photons. Such sources 
are more heavily biased, and their lower number d ensity 



yields larger and more spherical i onized regions (Zahn 
et al.||2007 [ [McQuinn et al.||2007b[ ). The same is true at 



low redshifts it feedback effects limit the role of low-mass 
sources. For instance, low-mass sources might suffer from 
thermal feedback onto t heir host galaxies during early 
stages of re ionization (Hult man Kramer et al.|2006 Iliev| 



et al. 2007). The effect on the patchy kti'L power spec 
trum of increasing the minimum source mass to 10^'^M(7 
(an increase of two orders of magnitude o yer the fidu 
cial Mmiii, discarding > 90% of sources, e.g., Trac fc Cen 



2007b [) is shown in the black solid and dashed lines of 



Figure [2] The impact of source bias on the shape of the 
patchy kSZ template is found to be small compared to 
our experimental uncertainty. 

Furthermore, recombinations in dense systems might 
effectively stall bubble growth early on in reionization, 
leading to a slo wly percolating web of similarly sized 
ionized regions (Furlanetto & Oh 2005). To estimate 
the effect of such a scenario on the kSZ power spectrum, 
we impose a relatively extreme maximum bubble scale 
of fmax = 5 Mpc/h. The result is shown in the blue 
solid line of Figure [2j Again, the changes are modest 
when compared to a model with the same duration but 
no limit on the bubble scale (dotted curve). 

We conclude that the shape of the kSZ power spec- 
trum due to patchy reionization is relatively robust to 
the effects of redshift translation, feedback, and recom- 
binations. Most importantly for this work, the amplitude 



is expected to be proportional to the duration of rcioniza- 
tion for uncorre latcd ionized regions and fixed Zmid (e.g., 
Gruzinov & Hu 1998). If the duration doubles in this 
Foisson process, there are twice as many independently 
moving ionized regions along the line of sight, which dou- 
bles the power. We find this to be a good approximation 
to the results from our simulations on the angular scales 
of interest (compare solid versus dotted lines in Figure 

Given the robustness of the shape of the patchy kSZ 
power spectrum to all these changes, we will assume a 
fixed template shape to derive our reionization results. 
As base template we choose a model with an efficiency 
of C = 20, making this model extend from Zx^=o.2 — 
11.0 to 2^^=0.99 — 7.8, however, as explained above the 
results are robust to variations in this choice. The base 
template used to compare to the SPT data was run at 
8 times higher spatial resolution than the simulations in 
the database. The resolution of the latter was chosen to 
lead to convergence at £ — 3000 to allow comparison to 
the data -derived constraint on the power on that scale, 
see S|5]2| 

3. DATA 

Here we describe the data used in this work, which 
serve to constrain either the ACDM model and primary 
CMB anisotropy, or the cosmic infrared background 
(GIB) and the thermal (tSZ) and kinetic (kSZ) Sunyaev- 
Zel'dovich effects. For the AGDM cosmological model, 
we include measurements of the baryon acoustic oscil- 
lation feature from SDSS and 2dFGRS (Percival et al. 
2010[ ), a measurement of the Hubble consta nt from the 
Hubble Space Telescope (Riess et al. 2011), an d mea- 
surements o f the prima ry GMH from W MA P 7 (Larson 



et al. 



20111. Most im- 



2011[ ) and SPT ( [Keisler et al 

portantly for this work, the CMB-|-hlAu-|-Ho preferred 
optical depth is r = 0.085 ± 0.014. The optical depth 
constraint is primarily due to the large-scale polariza- 
tion and polarization-temperature cross correlation data 



6 



Mm|n = M(Tvir=le4K), Zmid = 7.5, Az = 2 
Mmin = lelO Msol, Zmid = 7.5, Az = 2 
Zmid = 7.5, Az=4 
Zmid = 11.5, Az = 2 
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Fig. 2. — kSZ power spectra due to patchy rcionization for differ- 
ent Zmid, defined by Xe = 0.5 (dot-dashed versus dotted black 
lines), minimum star forming mass Mmin (dashed versus solid 
black lines), and maximum ionized bubble scale (blue solid 
line). The different Mmin cases are shown at a smaller Az to al- 
low plotting in the same panel. The blue solid line shows the kSZ 
power spectrum for a reionization model with Tmax — 

5Mpc/h. 

The red, dot -dash ed line shows our fiducial homogeneous kSZ 
template (see i ]4.f 



from WMAP7. 

The above datasets do not constrain the Sunyaev- 
Zel'dovich effects or CIB, which dominate the 
anisotropics at arcminute scales and SPT wave- 
lengths. For these small scale anisotropics, we include 
new multi-frequency data from the SPT (Rll) and the 
recent Planck /KFl C IB bandpowers at 217 GHz ( [Planck 
Collaboration 20111. The SPT bandpowers are based 
on observations at 95, 150, and 220 GHz of 800 deg^ of 
sky. From these maps, Rll estimate multi- frequency 
bandpowers covering angular scales 2000 < £ < 9400. 
The three-frequency SPT bandpowers are used to 
constrain the CIB, tSZ, and kSZ power spectra. The 



Planck /HFl data aid in constraining the CIB frequency 
dependence and its angular dependence on large scales. 

4. MODEL FITTING 

We use Monte Carlo Markov chain (MCMC) tech- 
niques to fit the data with a ACDM cosmological 
model including lensed primary CMB anisotropy, tSZ 
anisotropy, kSZ anisotropy, and foregrounds. The model 
and parameter fitting is described fully in Rll. Here we 
briefly describe the full model before focusing on the ele- 
ments shown by Rll to be important for the kSZ power 
measurement. 

We adopt the standard, six-parameter, spatially flat, 
ACDM cosmological model to predict the primary CMB 
temperature anisotropy. Rll consider model extensions 
including running of the spectral index, massive neutri- 
nos, and freedom in the number of neutrino species, find- 
ing no effect on the kSZ constraints. Rll also show that 
the systematic uncertainty in kSZ power due to uncer- 
tainty in the CMB lensing potential is much smaller than 
the derived kSZ uncertainties. 

We include seven free parameters for the small-scale 
temperature anisotropics from the tSZ and kSZ effects, 
and the CIB. We also include radio galaxies and galactic 



cirrus with strong priors as described by Rll. Two pa- 
rameters are the amplitudes of the tSZ and kSZ power 
spectra; the remaining five parameters describe the CIB 
mode l. The chosen tSZ template is the baseline model 
from Shaw et al. (2010). Rll show that the resulting 



kSZ power spectrum results are insensitive to the shape 
of the tSZ template. The kSZ template and CIB model 
are critical for the reionization constraints in this work, 
and will be described next. 

4.1. Kinetic Sunyaev-Zel'dovich effect 

We adopt two main templates for the kSZ power spec- 
trum. The first template includes only the contribution 
from patchy reionization (see fpl. The second and more 
realistic template includes roughly equal contributions 
from the reionization and post-reionization epochs. We 
henceforth refer to the post-reionization, homogeneously 
ionized, component as the "homogeneous kSZ" signal, 
and that from reionization as the "patchy kSZ" signal. 
We discuss the homogeneous kSZ model next. 

4.1.1. Homogeneous kinetic Sunyaev-Zel'dovich effect 

For the homogeneous kSZ signal, we adopt the cool- 
ing plus sta r-formation (CSF) model presented by Shaw 
et al. (2011 ). This model is constructed by calibratmg an 
analytic model for the homogeneous kSZ power spectrum 
with a hydrodynamical simulation including metal licity 
dependent radiative cooling and star-formation. [Shaw 



et al. (2011 ) measured the power spectrum of gas density 
fluctuations in the simulation over a range of redshifts, 
and used this to calculate the homogeneous kSZ power 
spectrum. The approximate scaling for this model with 
cosmological parameters is given by 



D 



homog 
3000 



1.9 fiK^ 
0.264 / 



h 
071 

-0.44 



1.7 

Vo.80/ 



4.5 



0.044 



096 



-0.19 



(10) 



This scaling (but not amplitude) is also a good approx- 
imation to the cosmological dependence of the patchy 
kSZ component. We vary the predicted post-reionization 

Djqqq"^ with Zond when interpreting the results in <|6| The 
CSF model kSZ component for Zend = 8 is shown m the 
red line in Figure [2| The model utilized here assumes he- 
lium is singly ionized between 3 < z < Zend and doubly 
ionized for z < 3. 

4.2. Dusty galaxies (CIB) 

The CIB is produced by thermal emission from dusty 
star-formin g galaxies (DSFGs) over a very broad range 



Marsden et al.[ [2009) 
Tom a tew molecules to 



in redshift (Lagache et al. 2005 
The dust grains, ranging m size 
0.1 mm, absorb light at wavelengths smaller than their 
size, and re-radiate it at longer wavelengths. Sufficient 
absorption occurs to account for roughly equal amounts 
of energy in the CIB and in the unproces sed starlight that 
makes up the o ptical/UV background (Dwek & Arendt 
1998; Fixsen et al. 1998). 

'I'he power spectrum of these DSFGs will have Poisson 
and clustered components. Both CIB components have 
more power than the kSZ power spectrum on the angu- 
lar scales and photon frequencies of interest (except at 
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95 GHz where the powers are comparable) . The Poisson 
component shows up as cx since the galaxies are 
small (effectively point sources) compared to the angular 
scales probed here. We model the clustered component 
by a power-law of the fo rm Dj^ oc i'^'^, a shape motivated 
by recent observations (Addison et al. 2012 Rll). We 
have explored allowing freedom in the angular shape via 
freeing the power-law exponent or adding a linear-theory 
template, and find the additional shape parameter in- 
creases the kSZ uncertainties by ^10% without changing 
the mean. 

The frequency dependence of the GIB is a key question 
for the kSZ measurement. Essentially, the GIB is well 
measured at > 220 GHz by SPT, Planck, and other ex- 
periments, but the signal-to-noise falls off towards longer 
wavelengths. The GIB power at 95 GHz is an extrapola- 
tion based on a fit to the modeled frequency dependence. 
Following previous work, we adopt the phenomenological 
modified black-body model: 



(U) 



V is the observing frequency, and (T) is the black-body 
spectrum for temperature T. T and /? are free parameters 
with uniform priors T G [5, 35 ivT] and P e [0,2]. Rll 
also e xamined the single-SED model presented by |Hall| 
et al. ( |2010) . The resulting kSZ power upper limits were 
similar when allowing for the possibility of correlation 
between the GIB and tSZ and 20% tighter when any 
correlation was assumed to be negligibly small. 

We assume the clustering and Poisson components 
have the same frequency scaling in this work. Implicitly 
this is equivalent to assuming that the s ource popu l ations 
and redshift distributions are similar. Hall et al. (20101 
argued that any difference in the frequency scaling should 
be small based on GIB simulations. As mentioned in Rll, 
the data have some preference for differing spectral in- 
dices when considering the 95 - 220 GHz alone, however 
this preference vanishes when Planck /HFI 353 GHz data 
is added. 

4.3. Correlation between the tSZ and CIB 

An important component of our model is the poten- 
tial for a spatial correlation between the DSFGs that 
produce the GIB and the groups and clusters that gen- 
erate the tSZ power spectrum. At ^ = 3000, most of 
the tSZ power is predicted to originate from clusters of 
mas s greater than 5 x 10^^ Mp) and redshift less than 
1.5 dShaw et al.||2009l iBattaglia et al.||2011l iTrac et al 



201ip. It the galaxies residing in these halos contribute 



a significant fraction of the GIB power at the SPT fre- 
quencies, then the spatial correlation between the tSZ 
and GIB signal will be non- negligible. Below 220 GHz, 
the tSZ effect produces temperature decrements; we thus 
expect an anti-correlation between the GIB and tSZ fiuc- 
tuations. Note that cross-spectra involving the 220 GHz 
band - the tSZ null frequency - may still contain power 
due to tSZ-GIB correlation as the GIB signal couples with 
the tSZ at the second frequency. 

Following Rll, we model the contribution of a tSZ-GIB 
correlation to the measured power spectrum as: 



Here D^^^,^^^^ is the power due to correlations, f^'^fj,. 
is the tSZ power spectrum at frequency i, and Dfl^. is 
the sum of the Poisson and clustered GIB components. 
^ is the correlation coefficient; we define ^3000 as the am- 
plitude of this correlation at £ = 3000. We consider two 
templates for the spatial dependence of the correlation, 
one fiat and one rising with £. The reasoning behind 
these templates is discussed later in this section. 

Rll measure the magnitude of the tSZ-GIB correla- 
tion by introducing ^ as a free parameter in their model 
fit to the SPT band powers, assuming an ^-independent 
template for Df^''^^^. They find C3000 = -0.18 ± 0.12, 
consistent with a fraction of the GIB emission bein g spa- 
tially c oincid ent with the tSZ signal. However, George] 



et al. (|2011|) extrapolate the 24 //m flux of galaxies in 



their (X-ray selected) group catalog to 2 mm and find a 
very small signal. Even allowing for uncertainties in this 
extrapolation, they conclude that the contamination of 
DSFGs to the tSZ signal is no more than a few percent. 
This would imply only a very weak spatial correlation be- 
tween these signals. A robust measurement of the total 
DSFG emission in groups and clusters at mm/sub-mm 
wavelengths would place tight observational constraints 

There currently exist no published theoretical esti- 
mates of the magnitude of the tSZ-GIB correlation. How- 
ever, we find that the publicly available simulations of 



Sehgal et al. (2010) predict a large anti-correlation. Tak- 
ing the cross-spectrum of their simulated GIB-only and 
tSZ-only maps at 148 GHz, we find ^3000 = —0.37. 

We make new predictions for the tSZ-GIB corre lation 
by combining the tSZ model of Shaw et al. ( 2010[) with 
the halo model based GIB calculations ot ISnang et aL] 



lates 

mass function of Tinker et al. 



([2012 1. In brief, the [Shang et al.| (|2012 ) model popu- 
ralos and sub-ha ios with USFGs utilizing the halo 

j|2p08) and the sub-halo 



mass function of Wetzel fe White ( f2tll0| . It is assumed 



that the relation between galaxy luminosity and sub-halo 
mass is lognormal, with free parameters for the overall 
normalization, characteristic mass, and redshift evolu- 
tion. These parameters are calibrated by comparing the 
model predictions with the recen t Planck measurement s 
of the GIB p ower spectrum ( Planck Gollaboration|20lT I 



Shang et al. (2012) explore model parameter constraints 
tor various scenarios in which the characteristic dust tem- 
perature of DSFGs, the redshift evolution of the subhalo 
mass - DSFG luminosity relation, or the shape of their 
SED are alternatively fixed or allowed to v ary. For each 
of the se scenarios (labeled "cases to 5" in Shang et al. 



(12) 



2012), we calculate the correlation coefficient, ^ 

We find that the predicted value of ^3000 varies between 
—0.02 > ^3000 > —0.34, well within the range found by 
Rll. ^ depends sensitively on the assumed redshift evo- 
lution of DSFG luminosity (or equivalently, the evolu- 
tion of their star- formation rate), which always increases 
towards higher redshift. For models in which the lumi- 
nosity evolves rapidly, a larger fraction of the GIB power 
is contributed by high-redshift (2 > 2) objects and thus 
the magnitude of ^ decreases. 

The ^-dependence of ^ also varies between our GIB 
models. In Figure[3 | we plot £f fo r each of the model cases 
presented by |Shang et al. (2012). The values of ^3000 are 
given in the caption. I'he shaded region represents the 
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Fig. 3. — The correlation coefficient £ calculated for the var- 
ious CIB model cases presented in Shang ct al. (2012). The 
correlation coefficients at I = 3000 are (trom top to bottom) 
{-0.34, -0.09, -0.07, -0.04, -0.03, -0.02}. The slope of the cor- 
relation with £ becomes steeper as |£| decreases. The magnitude 
of £ is most sensitive to the redshift evolution of the halo mass - 
DSFG luminosity relation; a more rapid increase in luminosity with 
redshift shifts the peak of the CIB to higher redshifts and there- 
fore reduces the strength of the tSZ-CIB correlation. The shaded 
regions represent the 1 and 2 a constraints on £ presented in Rll 
based on an ^-independent correlation template, over the range of 
angular scales probed by SPT. In this work we also explore the 
effect, on the derived kSZ power, of using an ^-dependent template 
with a free amplitude, shown in the solid line. 

1- and 2 a constraints on ^ obtained by Rll. We find a 
general trend that the correlation coefficient is a stronger 
function of £ for models in which the overall amplitude 
is closer to zero. This is straightforward to understand; 
models with a lower value of |^| are those for which high 
redshift galaxies contribute a larger fraction of the total 
CIB power. Most of the tSZ-CIB correlation then comes 
from higher redshift halos which principally contribute to 
the tSZ power at small scales. Therefore, the correlation 
coefficient increases with £. 

We will explore the dependence of constraints on the 
kinetic SZ power spectrum on the shape of , focusing on 
the model shown in the solid curve, which has an average 
correlation coefficient of 0.04, but letting the amplitude 
freely vary. Note that this is conservative in that models 
with steeper amplitude have correlations close to zero, 
however we will not impose a prior. 

5. RESULTS FOR THE OPTICAL DEPTH AND PATCHY 
KSZ POWER SPECTRUM 

In this section, we briefly review the optical depth mea- 
surement from WMAP7 in the context of reionization. 
We then present the SPT patchy kSZ power spectrum 
measurement, and discuss factors that can affect this re- 
sult. 

5.1. Optical Depth 

WMAP measures the optical depth to rei onization to 
be T = 0.088 ± 0.014 ( [Larson et aT][20lT| . Figure ^ 
provides intuition on how this constrains the epoch 
reionization. In the left panel, we show ionization his- 
tories for several models from our simulation database. 
The right panel plots th e large-scale pola rization signal 
as calculated by CAME ( Lewis et al.|2000 1 from the sim- 
ulated ionization histories along with the WMAP 7 band- 



powers. The near instantaneous model at ^mid = 10-5 
shown in the purple dotted curve, as well as the more 
extended green dot-dashed curve {Az ~ 3) and very ex- 
tended blue dashed curve {Az ~ 8) are good fits to the 
WMAP data, as shown in the right hand panel of Figure 
|4j The large scale polarization data do not allow one 
to discern between shorter and longer duration models. 
However, the red dot-dot-dashed and orange solid curves 
with Zmid — 7 and 12.5 are poor fits to the WMAP data. 
For the data used in this paper (CMB-f BAO+Hq), the 
optical depth constraint is r = 0.085 ± 0.014. This cor- 
responds to a ~ 0.25 shift to lower redshifts compared to 
using the WMAP-only optical depth constraint; we use 
r = 0.085 ± 0.014 for all reionization constraints in this 
work. 

5.2. kSZ power 

The SPT constraints on the kSZ power spectrum are 
summarized in Figure |5] and Table [T] The table contains 
95% confidence upper limits on the patchy kSZ power at 
£ — 3000 for different homogeneous kSZ and CIB models. 

The kSZ contributions from various redshifts have the 
same frequency dependence, and only small differences 
in angular dependence (as shown in Figure l2|. We find 
that the 2008 and 2009 SPT data presentecTby Rll do 
not yet significantly discriminate between the patchy and 
homogeneous kSZ template shapes. As a result, the mod- 
eled homogeneous kSZ power affects the level of patchy 
kSZ allowed by the data. With more homogeneous kSZ 
power, less patchy kSZ power is allowed. We note, how- 
ever, that the patchy kSZ likelihood curves in Figure |5] 
shift by relatively little when changing the homogeneous 
kSZ model by ~ 2fiK^. This is because the likelihood 
curves peak at or below zero kSZ power (—1.2 a with- 
out tSZ-CIB correlation), and a positivity prior is ap- 
plied. As will be discussed below, the SPT data prefer a 
non-negligible tSZ-CIB anti-correlation together with an 
increased kSZ amplitude. 

We consider two models for the homogeneous kSZ sig- 
nal. The first, zero homogeneous kSZ (dashed lines in 
Figure |5]), is the most conservative when determining 
upper limits on the patchy kSZ signal. However, it is 
physically unmotivated. The secon d case (solid lines) is 
the CSF model described in ^ 4.1.1 which is our best es- 
timate for the homogeneous kSZ signal. We scale the 
homogeneous kSZ signal according to the cosmological 
parameters at each step in the chain following Eqn. |10[ 
We tabulate the patchy kSZ upper limits with these two 
homogeneous kSZ models in the first and second rows of 
Tabled 

To constrain the epoch of reionization, we use kSZ 
power constraints based on a fixed template set to the 
sum of the CSF and patchy kSZ models in the fiducial 
cosmology. Constraints on the total kSZ power with this 
template are found in the third row of Table IT] The 
small error caused by our choice of a fixed ratio between 
patchy and homogeneous kSZ templates amounts to, at 
most, a fraction of the observed 15% difference between 
the kSZ upper limits for the CSF and patchy templates. 

A significant uncertainty in the current kSZ constraints 
is related to t he unknown tSZ-CIB correlation. As dis- 
cussed in |4.3[ the amount of correlation is highly uncer- 
tain; we therefore explore three cases for the correlation. 
For no correlations, shown in the blue lines of Figure [5] 
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Fig. 4. — Ionization history and integrated optical depth for different reionization models. The color-coding of model lines is the same 
between panels. Left parted' Evolu tion of the ionized fraction for each model. Right panel: EE power spectra for the models, together 
with WMAP7 ( [Larson et al.|[201l[ | band power constraint (black data point). The WMAP measurement of E-mode power in the band 
2 < £ < 7 does not restrict the duration of reionization. 

TABLE 1 

h.kSZ, ^ G [-1,1] kSZ Upper Limits 

CSF h.kSZ, ^ G [-1,1] 
h.kSZ, 5| G [-1,1] 
CSF h.kSZ, 5| G [-1,1] 
h.kSZ, ^=0 
CSF h.kSZ, ^=0 

Note. — 95% confidence upper limits on the kSZ power at 
I = 3000 for different assumptions about the homogeneous kSZ 
model and tSZ-CIB correlation, ^. The first and second rows have 
limits on the patchy kSZ power assuming either zero homogeneous 
kSZ power or the CSF homogeneous kSZ model. The third row 
shows constraints on the total kSZ power for a combined patchy + 
CSF homogeneous kSZ template. The combined template is used 
to infer constraints on the epoch of reionization. The first column 
assumes the most restrictive model for the correlations, § = 0, 
and therefore leads to the tightest kSZ limits. The second colu mn 
allows free correlations with the rising model shape (see j ]4.3| and 
the solid line in Figure |3| while the third column allows a free, 
^-independent tSZ-CIB correlation. The latter yields the weakest 
kSZ upper limits. 

contribution, and tlie tSZ-CIB correlation, f . The fig- 
ure assumes no liomogeneous kSZ contribution and an £- 
independent tSZ-CIB correlation; qualitatively the result 
is independent of these assumptions. It is evident that 
the kSZ power increases with the anti-correlation (more 
negative values of the correlation, ^). Anti-correlation 
has the opposite effect on the tSZ power. The decrease 
in tSZ power with increasing anti-correlation is shown 
in the color-coding of the independent samples, ranging 
from values of 4 fj,K^ for the tSZ with zero tSZ-CIB cor- 
relation, to small values of 2 fiK^ for a -0.5 correlation. 

Decreasing tSZ power with increasing anti-correlation 
contradicts our naive expectations. In single-frequency 
band powers (e.g., only 150 GHz), increasing tSZ-CIB 
anti-correlation would increase the allowed tSZ contribu- 
tion for any observed power level. The anti-correlation 
term is negative in this band, thereby allowing more tSZ 
power. However, this picture changes when we consider 
the SPT multi- frequency data. In Figure [7j we show rel- 
evant model components across the three deepest SPT 
frequency combinations (95 x 150 GHz, 150 x 150 GHz, 
and 150 x 220 GHz from left to right). We have nor- 
malized each model component to 2 fiK^ at 150 GHz and 
£ = 3000. Note that a tSZ-CIB anti-correlation would 



Fig. 5. — Constraints on the amplitude a.t i = 3000 of patchy 
kSZ fluctuations from the SPT data for different homogeneous kSZ 
and tSZ-CIB correlation models. The upper limits tighten with 
the addition of homogeneous kSZ power and with more restrictive 
assumptions about the tSZ-CIB correlation. 

and left column of Table IT] the inferred 95% confidence 
upper limits on the patcfy^ kSZ power a.t £ — 3000 are 
2.5 and 2.1 /iK^ with the zero and CSF homogeneous 
kSZ models respectively. Our second patchy kSZ con- 
straint (black lines in Figure [5] and center column of Ta- 
ble [l]) uses a risi ng t SZ-CIB correlation shape based on 
the modeling in j ]4.3[ This correlation shape is the solid 
hue in Figure [3] and corresponds to the ^3000 = —0.04 
Shang model case. However, the amplitude of the tSZ- 
CIB correlation is allowed to vary freely in the JVICIVIC. 
The third case has an ^-independent tSZ-CIB correlation 
of unknown amplitude. The data prefer a larger anti- 
correlation in this final case; the 95% confidence range is 
C3000 G [—0.28,0.03] for the rising correlation shape and 
■Caooo S [—0.43, —0.01] for the ^-independent correlation 
shape. As a result, the third case yields the most conser- 
vative constraints on the kSZ effect. The 95% confidence 
upper limits on the patchy kSZ power weaken to 5.5 and 
4.9 /xK^ with the zero and CSF homogeneous kSZ models 
respectively. 

Figure |6] shows the degeneracy between the patchy kSZ 
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Fig. 6. — Degeneracy between the patchy kSZ and tSZ-CIB cor- 
relation fit to the SPT data. The dots are de-correlated Markov 
chain samples, color coded by the tSZ power. The inner/outer con- 
tours are 68% and 95% regions, respectively. There is little room 
for patchy kSZ power with § = 0. Allowing ^ to be free, the patchy 
kSZ contribution can be large at the cost of small values of tSZ 
power. 

lead to a negative power rather than the positive sig- 
nal (shown in the blue dot-dashed line) in the plot. As 
shown, a tSZ-CIB anti-correlation reduces the power at 
150 X 150 GHz (which could be compensated by increas- 
ing tSZ power, shown in the black dashed line) but it 
also has ~50% larger effect at 150 x 220 GHz (where 
there is effectively zero tSZ power). Conversely, a tSZ- 
CIB anti-correlation reduces the power at 95 x 150 GHz 
by a similar amount to 150 x 150 GHz, whereas there is 
more tSZ power in the 95 GHz band. Adding kSZ power 
(shown as the combined template in the black solid line) 
more effectively compensates for the tSZ-CIB correlation 
term than adding tSZ power in the most sensitive com- 
binations of the three frequencies. A combination of the 
kSZ and tSZ does even better, as qualitatively shown by 
the red dotted line for 1.2Df^ - 0.2Df'^. 

6. CONSTRAINTS ON THE EPOCH OF REIONIZATION 

In this section, we present constraints on the duration, 
beginning, and end of reionization. We first interpret 
the kSZ power in terms of the duration of reionization. 
We then combine the kSZ constraint from SPT with the 
integrated optical depth constraint from WMAP7 to con- 
strain the evolution of the ionized fraction as function of 
redshift, Xe- Finally, we use the measurement of Xg to 
present limits on Zcnd and Zbog- 

We derive constraints on the history of reionization by 
comparing the experimental constraints on the total op- 
tical depth and the kSZ power to the predicted values 
for each model in our simulation database. Recall that 
we use the reionization simulations to predict the patchy 
kSZ contribution. We then add this patchy power to the 



power predicted by the CSF homogeneous kSZ model at 
the fiducial cosmology for the Zend of this simulation to 
estimate the total kSZ power. Uncertainties in the cos- 
mological parameters within the Monte Carlo Markov 
chains used in this paper propagate into a 12.5% uncer- 
tainty in the kSZ power spectrum amplitude according 
to Eq. 10 We therefore include a 12.5% Gaussian scatter 
in I?3Q(^Q when deriving reionization constraints. This is 
a good approximation as the kSZ power constraint from 
the current data is independent of the CMB-derived cos- 
mological parameters. 

We marginalize over models in our simula- 
tion database, weighting each model by the joint 
WMAP+SPT likelihood. We use a conservative prior 
of = 1 at z < 5 based oii the model-independ ent 
constraints of ( [Fan et al.||2006b[ [McGreer et al.||2011 ). 



6.1. Constraints on the duration of the epoch of 
reionization, Az^i 

Here we present constraints on the duration of reion- 
ization, defined as Azj-oi = %^=o.20 ~ %e=o.99- These 
constraints are primarily due to the SPT kSZ measure- 
ment; WMAP data alone do not discriminate between 
reionization models of varying durations. However, by 
constraining the integrated optical depth, the WMAP 
data indirectly affect the duration inferred from the SPT- 
derived kSZ constraint since models with given duration 
have a larger signal at higher redshift due to the higher 
mean density at higher redshift (see Figure 2|. 

In Figure [Sj we show constraints on Az f or the three 
tSZ-CIB correlation cases described in [5.2 Under the 
assumption of no tSZ-CIB correlation, we find a 95% 
upper limit on the duration of Az < 4.4. The rising 
correlation shape template leads to an intermediate limit 
of Az < 6.2. With the ^-independent correlation, we find 
the weakest 95% upper limit of Az < 7.9. Note from 
the left panel of Figure H that a shorter redshift interval 
follows the midpoint of ^he epoch of reionization than 
precedes it. As a result, only a fraction of the duration 
in Figure |8] can lie at lower redshift than the midpoint 
z - 10.25. 

6.2. Constraints on the evolution of the ionized fraction 

Xe 

In this section, we show that we can use the CMB data 
to constrain the history of reionization, Xe, by combining 
the constraints on the integrated optical depth and dura- 
tion of the epoch. To this end, we compute the posterior 
distributions of z{xe) for all values of Xe using our simu- 
lation grid given the constraints on t and the amplitude 

of the kSZ power spectrum, -Daooo'^^- 

As previously mentioned, the integrated optical depth 
constraint alone can not constrain the evolution of the 
ionized fraction as a function of time. Figure [9] shows 
the 68/95% confidence intervals in the Zcnd-Az plane. 
There is a perfect degeneracy in the large-scale WMAP 
polarization data between the duration and end of reion- 
ization; a low redshift for the end of reionization can be 
made consistent with the data by increasing the duration 
and the redshift of the starting point of reionization. The 
SPT constraint on the patchy kSZ contribution breaks 
this degeneracy by constraining the duration, as shown 
for the case of no tSZ-CIB correlations in the inner con- 
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Fig. 7. — Key model components in the three most important SPT frequency combinations. All model components have been normalized 
to equal Di = 2 ^K^ at ^ = 3200 in the 150x150 GHz spectrum. From left to right, the panels show the power spectra for the 95 X 150 GHz, 
150 X 150 GHz, 150 x 220 GHz bands. The SPT uncertainty on the average power across the multipole range £ £ [2600, 4600] is displayed by 
the black error bars (also centered at Di = 2 fiK'^ at £ = 3200). The tSZ-CIB cross power (blue dot-dashed line) increases slightly with 
frequency, as does a linear combination of tSZ and kSZ (red dotted line). As the tSZ (black dashed line) has the opposite frequency 
dependence, a tSZ-CIB anti-correlation allows for a larger kSZ amplitude (shown as the combined patchy+homogeneous template in the 
black solid line). The clustered component of the GIB is shown as well (green dot-dashed line); it has a much steeper frequency 
scaling. 
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Fig. 8. — Constraints on the duration of reionization for three 
assumptions about the tSZ-CIB correlation; free amplitude and 
scale-independent correlation (solid green line), free amplitude 
and rising correlation with amplitude consistent with modeling of 
low correlation simulations (dashed black line); and with cor- 
relation set to zero (dotted blue line). We have assumed the 
GSF homogeneous kSZ model for the post-reionization kSZ signal 
throughout. The WMAP large-scale polarization data enter by 
constraining the total electron scattering optical depth but are not 
able to constrain the duration of reionization. 

tours of Figure [9j 

To obtain constraints on the evolution of the ionized 
fraction with redshift, we integrate the posterior distri- 
butions at each Xf, to obtain the 68% and 95% hkehhood 
intervals for that ionization fraction to fall within a cer- 
tain redshift range. Again, we assume the CSF homoge- 
neous kSZ model for the post-reionization kSZ power. As 
indicated by Table [TJ the effect of different homogeneous 
kSZ power is modest given our tight upper limits. The 
constraints on Xe from the addition of SPT to WMAP 
are shown for the three tSZ-CIB correlation variants in 
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Fig. 9. — Gontour plot of z^^^ — Az plane showing 68% and 95% 
confidence regions. There is a degeneracy in the WMAP data be- 
tween the end of reionization and the duration (shaded contours) 
which is broken by the addition of SPT data (line contours). The 
tSZ-GIB correlation is assumed to be zero. 



Figure [To] The case with arbitrary, £- independent corre- 
lation (right panel, green lines) shows the weakest con- 
straints. The case with rising correlation shape (center 
panel, black lines) gives intermediate constraints. The 
case without correlations (left panel, blue lines) shows 
the tightest constraint on the ionized fraction, with the 
error on the timing dominated by the WMAP optical 
depth uncertainty. The bulk of the epoch of reionization 
is confined to a narrow redshift interval from z ~ 12.5 to 
z ~ 7 at 95% CL. 

We next consider the end of reionization which has 
been the focus of most current reionization observations. 
Recall that we define the end of reionization as Xe = 0.99. 
Constraints on Zcnd are shown in the left panel of Figure 
[TT] and tabulated in Ta ble [2] This is equivalent to taking 
a slice through Figure [Tofat Xe = 0.99. The SPT data 
enable the first CMB constraint on Zend- The SPT data 
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« = o 


e [-1,1] 


?e [-1,1] 


Az 


(68%) 
(95%) 


2.1 
4.4 


2.9 
6.2 


4.4 
7.9 


-^beg 


(68%) 
(95%) 


10.7 
12.1 


10.9 
12.3 


11.2 
13.1 


^end 


(68%) 
(95%) 


9.2 
7.2 


8.9 
6.4 


8.3 
5.8 



Fig. 10. — CMB constraints on the evolution of the ionized fraction for three cases of the foreground modeling. The SPT+WMAP 
68/95% confidence intervals are indicated by the thick/thin curves. The left panel shows the case without a tSZ-CIB correlation, the 
center panel shows the free ^£ case, the right panel shows the free fiat ^ case. 

disfavor a late end to reionization due to their prefer- 
ence for short durations. For ^ = 0, we find that hav- 
ing reionization end at z^nd < 7.2 is disfavored at 95% 
CL (inner blue contour, dotted line). The constraint 
weakens considerably when an arbitrary, ^-independent 
correlation is allowed (wider green contour, solid line). 
However, in this most conservative interpretation of our 
data, reionization still concluded at Zcnd > 5.8 at 95% 
CL. The rising-^ correlation shape yields an intermedi- 
ate constraint of Zend > 6.4 (cyan contour, dashed line). 
Note that while WMAP does not place a lower limit 
on the end of reionization, it trivially places an upper 
hmit (the WMAP contour hes directly beneath the SPT 
contours on this side) by limiting the integrated optical 
depth. Note also that the innermost contour is shifted 
slightly to the left compared to the cases with tSZ-CIB 
correlation. This is because the tighter SPT upper limit 
on the combined kSZ means that short models ending at 
lower redshift are preferred because they entail a smaller 
homogeneous kSZ contribution. In other words, there 
is a slight tension between the SPT requirement for a 
small kSZ amplitude and the WIVIAP requirement for a 
relatively large integrated opacity. 

Finally, we present constraints on when the first ion- 
izing sources turn on and begin reionizing the universe. 
Again, we have defined the beginning of reionization by 
Xe — 0.20. The likelihood function for Zbeg is shown in 
the right panel of Fig ure [TT| This is equivalent to taking 
a slice through Figure [T0| at x^, — 0.20. Wc find that the 
combination of WMAP and SPT data rule out an early 
onset of reionization at z > 12.1 at 95% confidence in 
the ^ = case. When allowing for a free, ^-independent 
tSZ-CIB correlation, the 95% confidence upper limit in- 
creases to Zbeg < 13.1. Again, the WMAP optical depth 
constraint leads only to a lower limit on the beginning of 
reionization, see above. Again the ^ = likelihood peaks 
at slightly lower redshift, for the same reason as in the 

■2^end CaSB. 

6.3. Forecasts for SPT full survey, Planck, and Herschel 

In the near future, we expect improved measurements 
of the kSZ power from the full SPT survey, optical 
depth (from Planck)^ and CIB (from Herschel and Her- 
schel /SW cross-correlation analyses). The SPT survey 



Note. — Reionization results for different tSZ-CIB correlation 
assumptions. We have defined volume-weighted ionized fractions 
of 20% and as the 99% as the beginning ^bog a^nd ^nd Zond of reion- 
ization, respectively, and the duration Az as the redshift interval 
between the two. We show 68% and 95% limits in each case as the 
likelihood surfaces are non-Gaussian. The first two rows (Az and 
•^bcg) show upper limits, while the third row (zcnd) shows lower 
limits. The first column assumes the most restrictive model for 
the correlations, § = 0, and therefore leads to the tightest reion- 
ization limits. The secon d col umn allows free correlations with the 
rising model shape (see i ]4.3| and the solid line in Figure [sjl while 
the third column allows a tree, ^-independent tSZ-ClB correlation. 
The latter leads to the weakest limits on the epoch of reionization. 

of 2,500 square degrees (three times the area used in this 
work) was completed in November 2011. The Planck sur- 
vey is ongoing and the first power spectrum results in- 
cluding the optical depth constraint from the large-scale 
E-mode polarization feature should be released in 2013. 
Herschel observations of the deepest 100 square degrees 
of the SPT survey will conclude in 2012, and should en- 
able a detailed study of the tSZ-CIB correlation as well 
as the CIB in general. 

We estimate the improvement in the kSZ power con- 
straint from the full SPT survey with 100 deg^ of Her- 
schel overlap by running a JVIonte Carlo IVIarkov chain 
with simulated bandpowers and uncertainties. We as- 
sume a 1% temperature calibration uncertainty and a 5% 
beam FWHIVI uncertainty in the SPT frequency bands; 
we assume the Herschel data is used to create a CIB tem- 
plate map which will be subtracted from the SPT bands 
to reduce the CIB contribution. Given the wide fre- 
quency range spanned, we introduce one new parameter 
to allow for decorrelation in the CIB between frequency 
bands. Instead of a single (3 in equation [TT] we assume 
a distribution of /3 in the galaxies of ^(P^p). We set 
a uniform prior of cr? G [0,0.35]; the upper edge is cho- 
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Fig. 11. — Likelihood functions for ^ond (^'e = 0.99, left panel) and ^bog of reionization {xe = 0.2, right panel). In both panels, WMAP 
data only requires that the end (beginning) of reionization will be after (before) the midpoint (yellow shading). The most conservative case 
with SPT data is shown by the green shaded region, solid line. This allows a tSZ-CIB correlation with an ^-independent shape. With the 
cyan region we show results when allowing correlations with the ^-dependent shape, dashed line. In the most restrictive case, we neglect 
correlations (blue shading, dotted line). We have assumed the CSF homogeneous kSZ model for the post-rcionization kSZ signal. 

sen such that the correlation between 150 and 220 GHz 
is at least 95%. For the full SPT survey without Her- 
schel, we find the kSZ constraint improves proportionally 
to the reduced SPT bandpower uncertainties; the uncer- 
tainties on the bandpowers and kSZ power are reduced 
by a factor of The combination of SPT and Her- 
schel leads to a factor of 6 improvement in the kSZ con- 
straints. Assuming the post-reionization kinetic SZ con- 
tribution is known, we find that the combination of SPT 
and Herschel data should be able to positively identify 
extended reionization at 95% confidence if Az > 2. Us- 
ing the Planck published sensitivity numbers leads to a 
forecasted constraint on the optical depth of At ~ 0.005 
( Zaldarriaga et al.|2008 ). This is almost three times bet- 
ter than th e current WMAP7- HSPT r =_a085 ± 0.014 

Figure 




constraint (Keisler et al. 2011) 



straint forecasts tor the end of reionization 



12 shows con- 
Our current 
.^-independent, tSZ-CIB 



constraints with an uncertain, 
correlation is reproduced for comparison in the green 
contour, solid line. We show three variations on the pre- 
dicted Zend likelihood function. The cyan contour (dotted 
line) shows constraints centered on the WMAP7-I-SPT 
optical depth if the combined future data set continues 
to set only a stricter upper limit on the kSZ power. The 
predicted Zend tightens significantly around Zcnd — 10.25 
and models that end at z < 8 would be ruled out at 95% 
confidence. At the same r value, a measured patchy 
kSZ amplitude of 2^K^, corresponding to a reionization 
duration of Az ~ 4 and in good agreement with the flat 
correlation at the preferred SPT value of ^ = —0.18 leads 
to the blue contour (dashed line). Reducing the central 
optical depth to the ~1 a lower bound of the WMAP 
measurement with zero tSZ-CIB correlation leads to the 
red contour (dot-dashed line). The latter two cases have 
similar effects on the Zend constraints. The predicted 
•Zend tightensaround Zond — 8 — 8.5 and models that end 
at z < 6.5 would be disfavored at > 95% confidence. 



7. DISCUSSION AND CONCLUSION 



Fig. 12. — Forecasted ^cnd likelihood functions for future 
ST'T+ Planck+ Herschel data. The blue region, dashed line 
shows the predicted likelihood function assuming zero tSZ-CIB cor- 
relation and the median WMAP7 optical depth. We explore the 
effects of shifting the best fit values for r to the WMAP7 lower 1 cr 
bounds in the red region, dot-dashed line, or alternatively of 
shifting the ^-independent tSZ-CIB correlation to the Rll prefer- 
ence of g = —0.18 in the cyan region, dotted line. The current 
conservative foreground modeling constraint is shown by the green 
region, solid line. 

In this work, we produce the first constraints on the 
evolution of the ionized fraction during the epoch of 
reionization using small scale CMB observations. To this 
end, we have presented a framework based on efficient 
reionization simulations to calculate joint constraints on 
Xe from measurements of the total optical depth and kSZ 
power. We have applied this method to new observations 
of the kSZ power from the SPT (R ll) and published op- 
tical depth results from WMAP7 ( [Komatsu et al.||2011 | 
to probe the beginning, end, and duration of the epoch 
of reionization. 

We find that the SPT kSZ measurement implies a short 
reionization duration. We show that the SPT kSZ con- 
straint is sensitive to modeling assumptions about the 
poorly known tSZ-CIB correlation. Assuming this cor- 
relation to be negligible leads to the strongest limits on 
the epoch of reionization duration. Conversely, assum- 
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ing no outside knowledge of the correlation leads to the 
most conservative results. The 95% upper limits on the 
duration are < 4.4 under the assumption of no cor- 
relations and Az < 7.9 for the most conservative as- 
sumptions. As a result, the SPT data combined with 
the WMAP optical depth constraint rule out reionization 
models that end very late or begin very early. With i- 
independent correlations, the epoch of reionization ends 
at z > 5.8 and begins at z < 13.1 at 95% CL. 

The end of reionization has been the subject of intense 
theoretical interest and observational effort. In Figure 
[T3j we contrast the most conservative CMB-derived con- 
straint on the neutral fraction 1 — Xe{z) (where Xe < 1) 
with other data. We choose to show l—Xe instead of Xg as 
this is more natural for the other datasets. Specifically, 
we are allowing an ^-independent tSZ-CIB correlation 
and are assuming the best-guess CSF homogeneous kSZ 
model for the post-reionization kSZ signal. The dark and 
light green contours are the 1 and 2 a likelihood contours 
respectively from the CMB data. Also shown in the up- 
per panel are ext ernal constraints on r eionization from 
the Ly-a fore st dMcGr ecr et al. 20111), Ly-a emitter s 
(LAEs) (Malhotra fc R hoads 20041 |Ouchi et al.|[20To|), 
a gamma ray burst GRB 050904 (Ibtani et al. 2006), 
and quasar proximity regions (Mesinger fc Haiman^2007j 



Mortlock et al. 2011). In the lower panel, we plot con- 
straints on the residual neutral fraction after reioniza- 
tion f rom the Ly-a forest ( Fan et aH2006b| jBecker et al. 



2007). The previously published data constrain the rel 
atively narrow interval of redshifts ~ 5 — 7. Most results 
are upper limits, meaning they are consistent with reion- 
ization having concluded much earlier. This is also true 
for the Ly-a forest points in the lower panel which cor 



respond to very lo w neutral fractions (Fan et al. 2006b 



Becker et al.|2007[ ). The SPT+WMAP7 data narrows in. 



on a previously unexplored region of the (z, 1 — a;e)-plane 
Two constraints are in mild tension with our results. 
These are the lower limits on 1 — Xe obtained from pos- 
sible detections of damping wing ab sorption in quasar 
proximity zones by Me singer fc Haim an (2004, 20 07[), 
and to a lesser exteniT Mortlo ck et al.| ( |201T| )7 |Boiton| 
et al. (2011). Addi tional proximity zone spec t ra, fur 



ther modeling (e.g. 



Bolton & Haehnelt 2007a Maselli 



et al.|[2007 Lidz et al.||2006 ), and constraints from other 



datasets should clarity the ionization state of the IGM 
at z = 6 — 7. 

In the future, it would be interesting to investigate 
whether there is tension with observations of the Ly-a 
forest after reionization. In particular, measurements 
of the mean transmitted flux through the Ly-a forest 
at z ~ 5 indicate that the ionizing sources emit only 
a few ionizing photons per hydrogen atom pe r Hubble 
time fMi ralda-Escude||2003[ fBolton fc Haehnelt||2007bj ). 
Since after accounting tor recombinations, a couple of 
photons per atom are required to complete reionization, 
these measurements suggest that reionization is a pro- 
longed process. 

The results presented here have implications for the 
proposed study of redshifted 21 cm emission from the 
hyperfine transition of neutral hydrogen. Observations 
of the hyperfine transition can potentially provide de- 
tailed three-dimensional information about the evolution 



new generation of telescopes will begin collecting data to 
detect this signa l. These eff orts include the MW A (Bow^ 
'man et al.ll2005l), LOFAR (Bar ker et aL] f20r2), GMRT 
{Ven et al. 2 009 J PAP ER (Pars ons et ^ 010) and SKA 
( Johnston et al.|2008 ). Improved constraints on the his- 
tory ot reionization can help these experiments select the 
optimal frequency range in which to focus their observa- 
tions. 

We also note that our results are consistent with the 
recent result based on the redshifted hydrogen hyperfine 
transition "global step" experiment EDGES,_whic h pro- 
duced a lower limit of Az > 0.06 ( Bowman fc Rogers] 
2010|. Our results are good news for upcoming experi- 



ments that hope to detect a redshifted 21 cm absorption 
feature sourced by substantial UV coupling p rior to X 



ray heating durin g the onset of reionization ( Pritchard 
fc Loeb|[20 10: Bow man et ar]|2008b|a[ [Burns et al.||^nnr 
Marker et al. 2012T The combination of SPT and 



and m orphology of the reionization process (e.g., Zaldar- 
riaga et al._^2004J . In the next few years, an aniibitious 



WMAP7 data limits the early stages of reionization to 
•^bcg ^ 13.1 at the 95% confidence level, meaning that an 
ionized fraction of 20% was reached later. While the hy- 
perfine structure absorption feature is expected at some- 
what lower ionization fractions, our result indicates that 
the redshifted 21cm line radiation from this epoch will 
appear at frequencies higher than those in models with 
higher Zbeg. There is less galactic synchrotron radiation 
and radio interference at these higher frequencies. 

With a wealth of new data on the horizon, this is an ex- 
citing period for CMB reionization constraints. The SPT 
survey of 6% of the sky is now complete. The kSZ con- 
straints from the full survey should be approximately \/3 
times better than those reported here; kSZ constraints 
will improve further with the deep SPTpol and ACTpol 
maps. The Herschel satellite is making unprecedented 
maps of the GIB, including a 100 deg^ overlap with the 
deepest 100 square degrees of the SPT survey. Further 
sub- mm surveys have sta rted observations or will do so 
in the near future (e.g., Wootten fc Thompson 2009 
Radford et al] [20071 [Holland et al.[|2006|) TEe cross- 
correlation analysis ot Herschel and SP'l' will provide 
information on the dominant foreground uncertainty for 
the kSZ power: tSZ-CIB correlations. Finally, the Planck 
survey is ongoing; the first power spectrum results are ex- 
pected in 2013. The Planck data should improve the op- 
tical depth constraint by a factor of 3. The combination 
of these data sets should be able to positively identify ex- 
tended reionization at 95% confidence if Az > 2, further 
illuminating the epoch of reionization. 
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Fig. 13. — CMB constraint on the redshift evolution of the mean neutral fraction for the (most conservative) free amplitude, ^-independent 
tSZ-CIB correlation case. The SPT+WMAP 68/95% confidence ranges are indicated by the thick/thin curves and the dark/light green 
shading. We also show other constraints on the neutral fraction based on quasar spectra (blue and violet constraints as well as black lower 
limits), a gamma ray burst (red upper limit), and hya emitters (green upper limits). 
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